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PREFACE 

Extensive  effort  has  been  instituted  both  in  the  U.S.  Coast 
Guard  government  efforts  and  in  the  commercial  community  for  the  further 
development  of  a  comprehensive  Loran-C  radio  navigation  system.  Vast 
areas  of  the  northern  hemisphere  are  now  covered  by  Loran-C  groundwave 
pulses  from  various  transmitters.  Receivers  for  position  location  and 
navigation  are  now  available  at  economically  reasonable  prices. 

Temporal  variations  in  the  propagation  time  of  the  pulses  due  to 
changes  in  the  electrical  boundary  condition  at  the  surface  of  the  earth, 
or  due  to  the  atmospheric  refractive  index  variations  resulting  from 
weather  systems,  interact  with  the  propagation  mechanism  for  the 
groundwave  and  result  in  degradation  of  the  navigation  accuracy  for  the 
navigator.  Temporal  propagation  effects  vary  with  geographic  locations, 
climate,  seasons,  and  perhaps  on  the  long  term  correlate  with  the  sun 
spot  cycle.  The  magnitude  of  these  effects  approach  the  order  of 
microseconds  in  anomalous  geographic  regions,  degrading  the  quality  of 
the  navigation  service  that  the  transmitting  system  provides. 

This  report  presents  the  results  of  a  theoretical  study  of  the 
possible  range  of  temporal  variations  in  the  100  kHz  groundwave 
propagation  mechanism.  This  report  includes  the  latest  state  of  the  art 
developments  in  LF  propagation  as  related  to  generalized  impedance 
parameters  and  parametric  evaluations  of  temporal  variations  related  to 
meteorological  parameters  of  the  earth's  atmosphere. 
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1. 


INTRODUCTION 


1.1  Loran-C  Description 

Loran-C  is  a  pulsed,  hyperbolic  navigation  system  operating  in 
the  90-110  kHz  band.  Ground  wave  range  is  typically  600  to  1400  nm  over 
seawater  and  considerably  less  over  land.  Predictable  accuracy  of  posi¬ 
tion  information  is  at  least  0.25  nm  (2-drms)  in  advertised  ground  wave 
coverage  areas  when  using  automatic  receivers  of  current  design.  The 
repeatable  accuracy  of  the  system  is  60  to  300  feet  (2-drms) ,  reference 
(1)  • 


Loran-C  has  been  adopted  by  the  U.S.  Coast  Guard  as  their  stan¬ 
dard  coastal  and  confluence  zone  aid  to  navigation.  By  1980-81,  the  en¬ 
tire  U.S.  Pacific  and  Atlantic  coasts  and  the  Gulf  of  Mexico  will  be  co¬ 
vered,  plus  the  Canadian  Pacific  Coast  and  the  Great  Lakes.  The  Canadian 
Coast  Guard  is  considering  a  further  Loran-C  expansion  over  the  Canadian 
Atlantic  coast  to  part  way  up  to  Labrador,  reference  (2). 

1.2  Accuracy  Requirements 

An  accuracy  of  0.25  nm  (=1500  feet)  is  normally  adequate  for 
navigation  on  the  open  seas  and  of  course  considerably  improved  accur¬ 
acies  are  achievable  as  the  distance  between  the  user  and  the  loran  chain 
stations  ^s  (because  of  the  higher  signal-to-noise  ratios) .  How¬ 

ever,  navigation  in  harbors  and  estuaries  requires  precise  accuracy  in 
order  to  navigate  safely  within  specific  channels,  especially  under  con¬ 
ditions  of  low  visibility  (i.e.  adverse  weather  conditions,  nightime)  and 
possibly  a  high  density  of  traffic. 

Therefore,  improved  accuracies  in  the  order  of  100  ft  (=100  ns) 
or  less  are  required  for  harbor  and  estuary  navigation.  To  accomplish 
this,  error  sources  must  be  identified  and  their  magnitudes  minimized. 
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1.3 


Sources  and  Magnitudes  of  Errors 


are: 


The  three  basic  sources  of  error  in  the  Loran-C  navigation  system 


•  Transmitting  System  Synchronization  errors 

•  Propagation  errors 

•  Receiver  measurement  errors 


Transmitting  system  errors  are  dependent  upon  the  stability  of 
the  frequency  standard,  transmitting  equipment,  timing  equipment,  output 
signal  phase  control  loops  and  the  system  area  monitor/synchronization 
techniques.  Propagation  errors  result  from  geophysical  and  meteorolog¬ 
ical  variations  along  the  signal  p>th ,  atmospheric  and  man-made  noise, 
and  electromagnetic  interference.  Receiver  measurement  errors  are  depen¬ 
dent  upon  receiving  equipment  performance,  antenna  environment,  relative 
motion,  and  system  geometry  (GDOP) . 


To  provide  increased  navigation  accuracy,  these  error  sources 
must  be  minimized.  In  doing  this,  the  three  modes  of  operation  of  a 
Loran-C  chain  should  be  considered.  These  are: 

1.  The  repeatable  modes,  i.e.,  ability  of  navi¬ 
gator  to  return  to  the  same  loran  coordinates 
again  and  again. 

2.  The  predictability  mode,  or  the  ability  of  the 
navigator  to  give  his  geographic  coordinates 
with  precision,  using  the  loran  coordinates. 

3.  The  relative  or  differential  loran  mode,  or 
the  ability  of  the  navigator  to  locate  posi¬ 
tion  relative  to  a  known  position  using  loran. 


Inaccuracies  introduced  into  the  system  due  to  transmitting  system  syn¬ 
chronization  errors  are  within  +  50  nanoseconds,  which  is  the  control 
tolerance  for  timing  adjustments  at  the  secondary  stations. 


1-2 


Receiver  measurement  errors  can  be  minimized  by  using  quality 
equipment,  proper  equipment  configuration  and  operation,  and  optimum 
crossing  angles  of  the  hyperbolic  lines-of-position.  A  theoretical  error 
budget  for  the  INTERNA V- LC2 0 4  hard  limiting  receiver  utilizing  100-second 
averaging  has  been  reported  by  GE-TEMPO.  reference  (3),  to  be  of  the 
order  of  10-17  nanoseconds.  Values  of  30-60  nanoseconds  error  for  a 
typical  commercial  receiver  using  5-second  averaging  are  also  reported, 
reference  (3). 

In  contrast  to  errors  of  the  order  of  say  50-100  nanoseconds  for 
the  transmitting  system  synchronization  and  the  receiver  measurements, 
propagation  errors  are  significantly  higher,  i.e.  400-3000  nanoseconds. 

Propagation  errors  occur  in  general  because  the  Loran-C  signal 
velocity  varies  with  position.  The  signal  is  slowed  down  by  the  physical 
and  electrical  properties  of  the  earth's/air  surface  interface.  These 
include  the  impedance  or  conductivity  of  the  ground,  the  roughness  or 
terrain  varations  of  the  surface,  the  refractive  index  of  the  atmosphere 
at  the  surface,  and  the  lapse  rate  or  (rate  of  change  of  refractive  index 
with  height  above  the  surface) . 

Spatial  (geophysical)  variations,  of  the  transmitted  signals  are 
primarily  influenced  by  the  inhomogeneities  of  the  surface  impedance  and 
the  variations  in  the  terrain.  Temporal  (meteorological)  effects  are 
produced  by  time  changes  in  these  spatial  features,  and  are  also  influ¬ 
enced  by  the  surface  refractive  index  and  the  lapse  rate  of  the  refrac¬ 
tive  index  of  the  earth's  atmosphere,  (which  are  known  to  change  diurnal- 
ly  and  with  changing  weather  conditions) . 

The  spatial  propagation  error  dominates  all  other  errors.  How¬ 
ever,  it  is  time  stable  and  its  effects  can  be  accommodated  through  grid 
calibration  programs.  Methods  involving  time  difference  measurements 
have  been  used  to  warp  the  Loran-C  hyperbolic  grid  lines  to  cancel  out 
these  fixed  spatial  dependent  errors. 

Temporal  related  propagation  effects,  although  of  smaller  magni¬ 
tude,  create  a  more  difficult  problem  for  maintaining  the  most  stable 
grid. 
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Thaoratical  atudiaa  aa  aarly  aa  l'»,'ft,  Hsfarancs  (4),  antlclpatad 
tha  axlatanca  of  tamporal  changaa  in  tha  phaaa  of  tha  Loran-C  signal  aa  a 
raault  of  natural  changaa  with  tlma  In  tha  propagation  machanlam  of  tha 
ground  wava.  Fortunataly  tha  ground  wava  axhtbita  a  ramarkabla  atability 
with  tlma.  Without  auch  atability,  Loran-C  praciaion  navigation  would 
not  ba  poaaibla.  Howavar ,  aa  tha  loran  tachnlcal  community  atrlvaa  for 
avar  graatar  accuracy  and  praciaion,  tha  amall  but  maaautabla  tamporal 
varlatlona  in  tha  ground  wava  propagation  machanlam  bacoma  mora  and  moi  a 
important.  With  a  total  allowabla  budgat  arror  in  tha  ordar  of  100  nano- 
aaconda,  tamporal  arror  a  ahould  ba  mlntmisad  to  about  10  nanoaaconda  for 
tha  combination  of  propagation,  transmittal,  and  tnoalvni  ayatam  arrora 
to  approach  100  nanoaaconda. 


OlVlKv'TtV*: 


2.1  Char acter iiation  of  Temporal  Propagation  Variation* 

The  primary  objective  of  thin  program  is  to  character  la*  tempo* al 
propagation  vat  iations  in  teima  of  secondary  phaae  cor  rect  l»>n*  aa  a  func¬ 
tion  of  distance  for  vaiytng  conditions  telated  to  meteor  ol»<g  lea  l  chang¬ 
es.  Ths  variable  pat  ametet  n  to  he  conn  l  det  ed  include  the  eutfav’e  i  eft  ac¬ 
tive  index,  the  vettical  lapse  factor,  and  the  electrical  chat actei 1st ics 
of  the  ground. 

Although  terrain  efTect  la  very  impoitant,  ite  effect  will  not  t»e 
considered  het  e  diiectly.  ilenetally,  1 1  t  e.jtt  l  at  1 1  les  in  tettain  increase 
the  effective  ground  surface  impedance,  x.  Theterote,  to  »xMiectly  in 
terpret  results  presented  hete,  the  electrical  chat actei  tat  tea  of  the 
ground  mentl>xred  alx've  should  he  taken  an  the  effective  am  face  imi'edance 
which  includes  tettain  effects. 


i.l  Assessment  of  Impact  on  lotan-0  Chain  Control 

A  secondary  ohlectlve  of  this  ptogtam  in  to  asses*  the  Impact  of 
using  knowledge  of  trmpotal  propagation  vat  iat  ivxts  to  improve  |,m  an-t* 
Chain  Control.  System  Atea  Mivtitora  t.'iAMs)  ate  lined  to  (nsttie  the  main 
tenance  of  the  emission  delay  at  the  sectindat  y  tiansmttteis  as  a  chain 
constant.  This  Imposes  problems  connected  with  sotting  out  system  t Imtng 
synchronlsat Ion  changes  ovet  the  baselines  and  propagation  changes  ovet 
the  path*  to  the  monitor.  The  system  timing  changes  .'vet  the  baselines 
have  been  obviated  by  the  use  of  the  t'esium  oactllatoi  clocks  Ixit  system 
changes  at  the  men i tot  still  terlect  both  system  timing  and  propagation 
changes  with  time.  Since  propagation  changes  ate  not  lineatly  distribut¬ 
ed  in  space  and  time,  the  temoval  of  any  gtons  system  changes  through  the 
dictates  of  the  atea  ithstltot  automat  teat  ly  tends  to  t  emove  the  nonlineat 
propaga t i on  temporal  changes  at  the  monitor.  However,  pi opagat l>xt  tern 
petal  change*  removed  at  the  tmxiltoi  geographic  l>H'*ti>Mt  will  not  lw»  ap¬ 
propriate  for  sixtte  other  point.  It  Is  therefore  important  that  propaga¬ 
tion  effects  in  the  time  domain  ttempotal  effects!  be  taken  into  account 
in  both  space  and  time  dimension*  as  the  chain  management  sttives  for 
ever  greater  Loran-C  accuracy  and  precisitxi. 


3.  ANALYSIS 


3.1  Definition  of  Variables  and  Causal  Mechanisms 


The  causes  of  ground  wave  temporal  changes  that  can  be  identified 
in  the  existent  theory  of  propagation  of  the  ground  wave  are  (not  neces¬ 
sarily  in  order  of  importance),  reference  (5): 


(1)  The  surface  impedance  of  the  ground  varia¬ 
tions  with  climate  and  weather; 

(2)  The  index  of  refraction  of  air  variations  at 
the  surface  of  the  ground  with  weather,  cli¬ 
mate  and  ground  elevation; 

(3)  The  variations  of  the  gradient  of  the  index 
of  refraction  at  the  surface  of  the  ground 
with  altitude  above  the  surface  of  ground 
(vertical  lapse  of  refractive  index);  and 

(4)  The  more  subtle  effects  of  terrain  roughness, 
elevation,  soil  consistency,  geologic  under- 
layment,  as  it  interacts  with  weather  and 
climate  to  produce  both  ground  impedance 
changes  and  atmospheric  changes  such  as  tem¬ 
perature  inversions. 


It  should  be  noted  that  temporal  effects  can  be  diurnal,  seasonal 
and  perhaps  even  exhibit  long  term  effects  correlated  with  the  sunspot 
cycle. 


3.2  Establishment  of  Variables  Base 


3.2.1  Surface  Index  of  Refraction  (na) 


The  index  of  refraction  of  air  at  the  surface  of  the  ground 
grossly  affects  the  primary  wave  propagation  time  and  is  usually  desig¬ 
nated  by  the  constant  value  n  ■  1.000338.  This  is  an  average  value 

a 
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given  in  reference  (6)  and  it  is  certainly  not  a  constant  either  geo¬ 
graphically  or  temporally.  Actually  na  may  vary  between  1.0002  and 
1.0004  or,  reference  (5),  between  200  and  400  N-units  where, 

N  -  (n-1)  106 

a 

Studies,  reference  (5),  show  that  the  secondary  phase  correction  is  sub¬ 
stantially  independent  of  the  surface  refractive  index.  However  there  is 
a  strong  correlation  between  the  surface  refractive  index  and  the  verti¬ 
cal  lapse  factor,  «,  and  a  strongly  affects  the  secondary  phase  cor¬ 
rection. 

3.2.2  Vertical  Lapse  Factor  (a) 

Changes  in  surface  temperature  and  indeed  surface  N-f actor  have 
been  correlated  strongly  with  the  vertical  lapse  factor.  Since  the  ver¬ 
tical  lapse  factor,  which  shall  be  designated  as  a,  enters  into  the 
ground  wave  theory  as  a  parameter,  calculations  of  the  ground  wave  param¬ 
etric  in  a  can  be  used  to  estimate  the  effects  of  climate  and  weather 
that  are  observed  at  the  surface  of  the  ground.  Of  course,  if  informa¬ 
tion  on  N-units  aloft  are  available,  these  also  can  be  used  to  deduce  the 
a  factor . 

The  temporal  variations  associated  with  weather  changes  at  the 
surface  of  the  ground  have  been  clearly  identified  in  a  number  of  papers. 
References  (7,  8,  9,  10  11  and  12).  Although  strongly  correlated  with 
the  surface  value  of  index  of  refraction  n(h),  h  -  O,  these  variations 
can  only  be  explained  by  the  atmospheric  vertical  lapse  variations  in  the 
index  of  refraction,  n(h),  where  h  is  the  altitude  above  the  surface. 

For  Loran-C,  a  may  vary  from  about  0.6  to  1.2  and  the  atmos¬ 
phere  n (h)  is  not  necessarily  exponential. 
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3.2.3  Ground  Electrical  Characteristics  (x) 

The  secondary  phase  correction  is  also  a  function  of  the  elec¬ 
trical  characteristics  of  the  ground.  These  ground  characteristics  can 
be  expressed  as  a  complex  number  x  which  is  a  normalized  surface  impe¬ 
dance.  If  an  effective  value  of  x  is  used  or  if  the  ground  is  homoge¬ 
neous,  the  classical  ground  wave  theory  permits  calculation  of  the  secon¬ 
dary  phase  correction.  Typical  values  of  x  are  given  in  Table  1  and 
Table  2,  reference  (5). 

The  values  of  amplitude  |x|  and  phase,  Arg(x),  given  in  Table 
1  define  the  region  of  the  complex  plane  that  has  been  found  through  ob 
servations  to  be  most  appropriate  to  Loran-C.  These  values  give  the  part  of 
the  complex  x-plane  that  has  been  used  in  practice  to  represent  the  ground 
surface  impedance  throughout  the  world  for  various  Loran-C  chains.  These 
values  cover  a  range  of  effective  values  that  are  physically  realizable.  In 
contrast,  ground  effective  conductivity  (o)  values  given  in  Table  2  may  or 
may  not  be  physically  realizable.  The  most  common  value,  0  -  5  mhos/m,  is 
valid  over  sea  water  and  corresponds  to  an  impedance  of  .001055  exp  (j 
.78535).  On  the  other  hand,  the  values  near  0  -  0.0001  mhos/m  are  not 
realized  at  Loran-C  frequencies  because  waves  penetrating  the  ground  as  the 
conductivity  decreases  are  reflected  at  various  levels  and  the  phase  of  the 
complex  impedance  (x)  moves  away  from  ff/4  (for  sea  water  )  to  typically 
inductive  values  of  impedance  near  1  radian. 

3.3  Secondary  Phase  Correction  Calculations 

The  phase  of  a  Loran-C  ground  wave  signal  is  usually  expressed  as 


TABLE  1 


RELATIONSHIP  BETWEEN  IMPEDANCE 
MODULUS,  j  x  |  ,  and  ARGUMENT,  Arg(x),  h  ■  0 


M 

— 

Arg  (x) 

1*1 

Arg  ( 

001 

.7854 

.04182 

.8198 

0010549 

.78535 

.043038 

.8268 

002 

.  7845 

.044305 

.8339 

003 

.7837 

.045 

.8377 

004 

.78  30 

.050 

.8635 

005 

.7823 

.055 

.8885 

006 

.7816 

.060 

.9150 

007 

.7809 

.08 

1.036 

008 

.7802 

009 

.7795 

010 

.7788 

.7781 

.012 

.7774 

013 

.7767 

.7760 

015 

.7753 

016 

.7745 

.7738 

018 

.7731 

Wm 

.7724 

1L1 

.7717 

021 

.7710 

.022 

.7702 

.023 

.7694 

024 

.7687 

.025 

.7679 

.025473 

.7672 

026979 

.7658 

.02719 

.7658 

027648 

.7658 

027793 

.7657 

028 

.7655 

029 

.7663 

030 

.7678 

.031 

.7700 

032 

.7730 

033 

.7762 

034 

.7800 

.035 

.7845 

.036 

.7890 

.037 

.7945 

.038 

.7995 

.039 

.8055 

.040 

.8115 

.04008 

.8116 

.041206 

.8163 

TABLE  2 


RELATIONSHIP  BETWEEN  CONDUCTIVITY  (O)  AND  HOMOGENEOUS 
GROUND  IMPEDANCE  WITHOUT  IRREGULARITIES  OR 
NON-HOMOGENEITIES 


o,  mhos/m 


.0001 
.0002 
.0003 
.0004 
.0005 
.0006 
.0007 
.0008 
.0009 
.0010 
.0020 
.0030 
.0040 
.0050 
.0060 
.0070 
.0080 
.0090 
.0100 
.0200 
.0300 
.0400 
.0500 
.0600 
.0700 
.0800 
.0900 
.1000 
.2000 
.3000 
.4000 
.5000 
.6000 
.7000 
.  3000 
.9000 


.11657 

.10470 

.09581 

.08882 

.08316 

.07846 

.07447 

.05273 

.04307 

.03730 

.03337 

.03046 

.02820 

.02638 

.02487 

.02359 

.01668 

.01362 

.01180 

.01055 

.00963 

.00892 

.00834 

.00786 

.00746 

.00528 

.00431 

.00373 

.00334 

.00305 

.00282 

.00264 

.00249 

.00236 

.00167 

.00136 

.00118 

.00106 

.00106 


tog  (*)ifk 


.42082 

.57579 

.64109 

.67588 

.69729 


where  a)  -  2irfj  (f  being  the  signal  frequency  in  Berts) 
na  -  refractive  index  of  air  at  ground  level 
c  •  velocity  of  light  in  vacuum 
d  "  transmitter  to  receiver  separation  distance 

<frc  -  secondary  phase  correction,  caused  by  finite  impedance 

and  refractive  index  discontinuities  across  the  boundary 
along  which  the  signal  propagates. 


Secondary  phase  correction  can  be  expressed  as: 


.  «/*  ,  L  „  1/3  2/3 

■  (alt.)  a  t 


o  \a 


-  radians 


or 


($</<*» 


10  Microseconds 


where 


a  ■  earth  radius 
k^  -  wave  number  ■ 


(1) 


is  primarily  influenced  by  surface  impedance  and  is  derived 
from  the  boundary  condition  at  the  surface  of  the  ground  as 
the  zero' th  root  of  a  Riccati  differential  equation  and  at  100 
kHz  is  typically  “0.42 


is  associated  with  the  lapse  rate  of  the  refractive  index. 
The  formula  relating  a  to  the  change  of  n  with  height  is 


a  -  1  *  1  5*  n  (a  »  h) 

nm  d  (a  h) 

cl 


(2) 


Equation  (2)  has  an  alternate  representation  in  terms  of  measur¬ 
able  atmospheric  variables. 


a 


N 

a 


) 


n  -  (n.  -  i)  io6 

a 

Nb  *  <nh  “  l>  10 

hfe  is  altitude  above  the  surface 


(3) 


where 


The  N-unita  at  any  height  can  be  written! 


N  -  (n-  1)  106  -  21 


**[»♦ 


together  with  the  hydrostatic  equation: 


where 


I  +  f  *  -10^ 


h  is  in  meters 

g  »  the  gravitational  constant  (cm  sec-2) , 
p  •  density  of  air  (grams  cnf^) , 

p  ■  atmospheric  pressure  £ millibar sj  (kg  cnT^sec-2) , 

T  ■  temperature  |°Kelvinj  , 

e  »  partial  pressure  of  water  vapor  ^millibarsj,  (kg  cm_1sec"2) . 

The  10  is  necessary  to  make  the  equation  dimensionally  correct  and  dimen¬ 
sionally  relative  to  equation  4. 


Then, 


dN  .  / 776  x  1.268  77.6  f 
dh  \  T  t2  [ 


9620e 


77.6  I  .  481ol  de 

T  [  “  T  J  dh 


Si 


where  a/n  “  a  and  a  *  6.378(10®)  meters,  (Reference  15.) 
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a  is  evaluated  at  various  temperatures  using  a  unit  height  of 
100  meters  and  P  equal  to  1013.25  mb.  Thus, 


at 


-15°C  a  •  .756  -  ^ .0753  +  .00277eJ  ^  +  .358  ^ 


at  0°C 


4  /  dT  de 

o  -  .770  -  <.0673  +  .  00234e>  ~  +  .301  ^ 


dh 


dh 


at  27  C  a  -  .791  -  <.0557  +  .00176e 


-f 


|  dT  t  ...  de 

dh  *264  dh 


at  35°C 


a  -  .796  -  | 


.0529  +  . 00163e  ^  ^  +  .251  ^ 


Also  for  homogeneous  atmosphere*,  de/dh  -  -  mg/RT  e  (Reference  16) 


«  -15°c  t 


at  0°C  - 

an 


at  27°C 

an 


« 35°c  f  ■ 


-. 0132e/100  m 

- . 0125e/100  m 

-.0113e/100  m 

-.0111e/100  m 


The  temperature  lapse  rate  wsuld  be  expected  to  be  between  the 
adiabatic  lapse  rate  of  -.98°C/100  m  and  the  dew  point  lapse  of  -  17°C/ 

100m.  Reference  16  gives  a  derivation  of  the  temperature  lapse  rate  in  a  homo¬ 
geneous  atmosphere  and  indicates  that  for  a  mixing  ratio  of  10  mille  which  cor¬ 
responds  to  16  mb  of  vapor  pressure  the  temperature  lapse  rate  would  be  one-half 
the  dry  adiabatic  value.  Therefore,  for  purposes  of  evaluating  a  values,  two 
average  temperature  lapse  rates  have  been  assumed  that  are:  -.98°C/100  m  and 
-.49°C/100  m. 


•Meteorological  explanation  of  homogeneous  atmosphere. 
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Over  land  paths  the  diurnal  temperature  will  vary  considerably 
du«  to  the  solar  heating  and  nocturnal  cooling  of  the  earth's  surface. 

The  extreaes  of  the  tesperature  variations  are  therefore  assumed  to  reach 
from  -4°C/100  a  at  noon  to  +4°c/100  a  at  midnight.  These  extreaes 
represent  a  8.16  or  4.08  multiplicative  factor  from  the  -.49  or  -.98°c/ 
100  a  lapse  rates.  The  water  vapor  lapse  rate  will  change  proportionally 
with  the  temperature  lapse  rate.  Therefore  the  saae  multiplicative 
factors  are  applied  to  the  de/dh  lapse  rate  as  are  applied  to  th  dT/dh 
lapse  rate  in  an  attempt  to  derive  a  values  for  winter  and  summer  con- 
ditons  from  the  formulas  given  in  the  previous  section. 

for  Winter  Conditions 


Moderate  Winter  0°C,  50%  humidity,  i 

.e.  e  » 

3  mb,^  -  -.0375 
an 

rab/100 

dT/dh 

de/dh 

Average  Conditions  -.98°C/100  m 

and 

-.0375  mb/100  m 

(Homogeneous  Atmosphere)  -.49°C/100  m 

and 

-.0375  mb/100  m 

Extreme  for  Noon  -  4°C/100  m 

and 

-.153  mb/100  m 

-  4°C/100  m 

and 

-.306  mb/100  m 

Extreme  for  Midnight  +  4°C/100  m 

and 

+.153  mb/100  m 

+  4°C/100  m 

and 

+.306  mb/100  m 

Extreme  Winter  -  15"  C,  50%  Humidity 

i.e.  e  ■ 

1  mb,^  -.0132 

rab/100 

dT/dh 

de/dh 

Average  Conditions  -.98°c/100  m 

and 

-.0132  rab/100  m 

(Homogeneous  Atmosphere)  -.49°c/100  m 

and 

-.0132  mb/100  m 

Extreme  for  Noon  -4°C/100  m 

and 

-.0539  mb/100  m 

-4°C/100  m 

and 

-.108  mb//100  m 

Extreme  for  Midnight  +4°C/100  m 

and 

+.0539  mb/100  m 

+4°C/100  m 

and 

+.108  mb/100  m 
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For  Summer  Conditlona 


Moderate  Sumer  27°C,  50%  Humidity,  i.e.  e  - 

17.5 

mb,|£  -  -.198  mb/100 

dT/dh 

de/dh 

Average  Conditions 

98°C/100  m 

and 

-  .198  mb/100  m 

(Homogeneous  Atmosphere) 

- . 49°C/100  m 

and 

-  .198  mb/100  m 

Extreme  for  Noon 

-4°C/100  m 

and 

-  .808  mb/100  m 

-4OC/100  m 

and 

-1.62  mb/100  m 

Extreme'  for  Midnight 

+4°C/100  m 

and 

+  .808  mb/100  m 

+4°C/100  m 

and 

+1.62  mb/100  m 

Extreme  Sumer  35°C,  High  Humidity,  and  e  • 

42  mb 

de 

■  -.466mb/100  m 

dT/dh 

de/dh 

Average  Conditions 

- . 98°C/100  m 

and 

-.466  mb/100  m 

(Homogeneous  Atmosphere) 

49°C/100  m 

and 

-.466  mb/100  m 

Extreme  Noon 

-4°C/100  m 

and 

-1.90  mb/100  m 

-4°C/100  m 

and 

-3.80  mb/100  m 

Extreme  Midnight 

+4°C/100  m 

and 

+1.90  mb/100  m 

+4°C/100  m 

and 

+3.80  mb/100  m 
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Calculated  a  values  are  given  below  in  Table  3.  These  can  be 
interpreted  as  follows: 

Winter  conditions  yield  inverse  correlation  between  propagation 
velocity  and  temperature  as  discussed  extensively  in  Reference  8.  That 
is,  the  lower  the  temperature  the  faster  the  signal  travels  or  the  smal¬ 
ler  the  phase  delay  for  a  given  path.  During  summer  conditions  a  nega¬ 
tive  term  associated  with  dT/dh  and  a  positive  term  associated  with  de/dh 
tend  to  cancel  each  other.  During  average  or  normal  summer  conditions 
these  factors  tend  to  offset  keeping  the  alpha  value  between  approx¬ 
imately  .75  and  .85.  Also,  under  extreme  summer  conditions  a  correlation 
with  temperature  can  average  all  the  way  from  the  winter  time  inverse 
correlation  to  a  high  temperature  summer  time  direct  correlation  shown 
for  +35°C  and  dT/dh  -  -.49°C/100  m  in  Table  3. 

TABLE  3 


a  Values  Calculated  for  Various  Conditions 
Listed  Above 


dT/dh  -  -.49  dT/dh  »  -.98 

Moderate  Winter  Conditions  are  T  *  0°C 


Average 

a 

=*  .795 

a 

• 

.832 

Extreme 

Noon 

a 

=■  .975 

a 

m 

1.02 

Extreme 

Midnight 

a 

-  .565 

a 

9 

.519 

Extreme 

Winter  Conditions 

are 

T 

9 

-15°C 

Average 

a 

=*  .790 

a 

M 

.828 

Extreme 

Noon 

a 

»  1.03 

a 

9 

1.05 

Extreme 

Midnight 

a 

=■  .482 

a 

9 

.463 

Moderate 

Summer  Conditions 

are  T 

'  a 

■  27°C 

Average 

a 

=*  .781 

a 

S» 

.823 

Extreme 

Noon 

a 

=*  .709 

a 

9 

.924 

Extreme 

Midnight 

a 

*  .873 

a 

9 

.658 

Extreme 

Summer  Conditions 

are 

T 

9 

35°C 

Average 

a 

-  .738 

a 

9 

.798 

Extreme 

Noon 

a 

=>  .328 

a 

9 

.804 

Extreme 

Midnight 

a 

-  1.26 

a 

9 

.787 
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Finally,  these  results  suggest  that  Loran-C  temporal  variations 
both  in  the  summer  and  in  the  winter  should  be  directly  predictable  from 
surface  measurements.  The  temperature  and  humidity  lapse  rates  dT/dh  and 
de/dh  will  be  determined  by  the  surface  values  of  the  temperature  and 
humidity  along  with  the  history  of  their  variations  as  observed  and 
recorded  in  the  immediate  past.  An  experimental  program  to  prove  this 
should  be  initiated. 

Another  example  derived  from  Reference  16  is  shown  in  Figure  1. 

The  information  in  Figure  1  and  in  Table  3  tends  to  show  the  same 
results.  Figure  1  shows  that  when  de/dh  in  mb/100  m  is  equal  to  .4  x  dT/dh 
in  °C/100  m,  the  changes  in  the  lapse  rate  dN/dh  will  tend  to  cancel. 

For  winter  conditions  where  de/dh  -  -.0132  mb/100  m  to  -.0375  mb/100  m 
and  dT/dh  •  -,49°C/100  m  to  -.98°C/100  m,  the  condition  for  equi¬ 
librium  is  not  met,  whereas  in  the  summer  when  de/dh  -  -.198  mb/100  m  to 
-.466  mb/100  m  and  the  dT/dh  «  -.49°C/100  m  to  ~.98°C/100  m,  the  con¬ 
dition  for  equilibrium  is  met  and  actually  reversed. 

All  of  this  material  then  tends  to  support  theoretical  predic¬ 
tions  of  Loran-C  temporal  effects.  It  supports  the  observed  large  winter 
temporal  variations  and  small  summer  temporal  variations  that  have  been 
consistently  observed  on  the  U.S.  East  Coast  path  between  Dana,  Indiana 
and  Carolina  Deach,  N.C.  An  experiment  utilizing  Loran-C  measurements 
and  surface  meteor  logical  measurements  over  paths  subject  to  extreme  sum¬ 
mer  and  winter  conditions  should  establish  this  relationship  through  data 
analysis  techniques. 
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A  N  Per  100  ■  A  N  Per  100  ■ 


4.  ANALYTICAL  RESULTS 


4.1  Atmospheric  Temporal  Effects 

Atmospheric  temporal  variations  are  correlated  with  surface  re¬ 
fractive  index  and  the  vertical  lapse  rate  of  the  surface  refractive  in¬ 
dex.  The  phase  temporal  variations  that  are  caused  by  change  in  n  , 
the  surface  refractive  index,  are  caused  by  changes  in  the  primary  wave. 
The  primary  wave  phase,  <t>,  is, 

<0  «  7  n  d,  (6) 

c  a 


and  the  phase  changes  associated  with  the  surface  refractive  index  is 
simply 


A  na*d, 

which  quantity  is  in  microseconds  if  d  is  in  kilometers.  Therefore,  a 
change  from 


na  -  1.0002  to  1.0004, 

(the  maximum  na  variation)  over  a  propagation  path  with  a  length  of 
1000  km  would  yield  a  phase  change  of  0.2  microseconds.  Normally  the 
change  in  average  value  (n  -  1.000338)  over  a  day  or  during  a  wea- 
ther  change  would  be  only  40  N-units  (N-units  =  (n  -1)  10  ) . 

Therefore,  the  maximum  diurnal  or  weather  change  effect  that  could  be 
attributed  to  changes  in  n  for  a  1000  km  propagation  path  would  be 

a 

0.04  microseconds,  or  40  nanoseconds.  This  amount  of  variation  may  be 
nearly  typical  of  the  variations  observed  over  a  1000  km  propagation  path 
during  warm  weather.  However,  in  cold  weather,  the  observed  variations 
are  typically  an  order  of  magnitude  or  more  greater  than  this  value,  re¬ 
ference  (8) .  As  a  result  the  overall  winter  weather  variations  observed 
over  an  approximately  1000  km  propagation  path  have  been  attributed  to 
changes  in  the  vertical  lapse  factor. 
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4.1.1  Graphs  Depicting  Atmospheric  Temporal  Effects 


Graphs  of  secondary  phase  correction  (tc)  vs  distance,  param¬ 
etric  in  vertical  lapse  factor  (a),  for  select  values  of  ground  surface 
impedance  x  (or  conductivity  0)  are  given  in  Figures  2  through  7. 

Since  there  is  a  negligible  effect  due  to  changes  in  the  surface  refrac 
tive  index  (n^j  except  in  so  far  as  na  is  correlated  with  a,  these 
curves  hold,  for  all  values  of  na  from  1.0002  to  1.0004.  a  values 
range  from  0.65  to  1.20  for  each  graph. 

Figures  2  and  3  are  very  similar.  The  former  is  based  on  a  sea¬ 
water  ground  impedance  x  *  .0011  (o  "4.6  mhos/m)  taken  from  Table  1 
in  Section  3  and  the  latter  is  based  on  seawater  ground  impedance  |x|  - 
.001055  (o  ■  5.0  mhos/m)  taken  from  Table  2  in  Section  3.  Although  not 
perceptable  on  the  graphs,  a  change  in  the  secondary  phase  correction  of  the 
order  of  nanoseconds  is  produced.  Furthermore ,  a  comparison  of  the  values  of 
the  secondary  phase  correction  using  seawater  ground  impedance  with  the  values 
in  Figures  4  through  7  deraomstrates  the  strong  dependence  of  the  secondary 
phase  correction  on  the  ground  surface  impedance.  Values  of  t^  at  1600  km 
over  seawater  (Figure  2  and  3)  are  of  the  order  of  2.5  to  4.0  us  whereas  for 
poor  earth  (Figure  7)  they  are  13  to  17  us. 

These  graphs  clearly  show  the  lmpot tance  of  a  in  temporal  prop¬ 
agation  changes.  The  change  in  t^  between  the  extreme  values  of  a  is 
>^1.5us  for  seawater  (Figure  21  and  «^4.5  us  for  poor  earth  (Figure  7) 
at  a  distance  of  1600  km.  It  should  be  pointed  out  that  the  phase  cor¬ 
rection  change  A( tc)  (g  proportional  to  a  A(c^/3)  and  a  Ad, where  d 
is  distance.  Since  A(a2/3)  is  smaller  than  Ad  the  effect  of  A(a2/3)  on 
A(tc)  is  more  evident  at  greater  distances.  An  expanded  graph  is  given  in 
Figure  0  to  more  clearly  illustrate  the  effects  of  a  at  short  distances. 

This  shows  that  At  is  about  O.lus  at  200  km  for  a  change  in  a  from  0.65 

V 

to  1.0  over  average  ground. 


4.1.2  Relationship  of  a  Values  With  Observations. 

The  seasonal,  diurnal  and  other  temporal  variations  observed  over 
the  1000  km  propagation  path  between  Dana,  Indiana  and  Carolina  Beach, 
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Figure  2  S>'ondary  phase  correction,  t  ,  microseconds,  as  a  function  of  distance,  d,  from  the  transmitter 
for  various  atmospheric  vertical  lapse  factors,  a,  and  a  ground  surface  impedance  for  sea  water 
x  =  .0011  exp  (j,7854)  or  a  conductivity  G  =  4.6  mhos/m  (Table  2). 
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Secondary  phase  correction,  tc,  microseconds,  as  a  function  of  distance,  d,  from  the  transmitter 
for  various  atmospheric  v«_rt  1  al  lapse  factors,  d  and  a  ground  surface  impedance,  x  =  .03337  exp 
< 3.77560),  and  ground  conductivity  0  =  .005  mhos/m  (Table  2). 
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Figure  8  tc  vs  Distance  for  Vertical  Lapse  Factor  (a) 

of  1.0  and  0.65,  for  Ground  Impedance  x  =  .033  exp  (J. 77620) 


reference  (3),  can  be  explained  by  an  a  value  of  1  +.02  during  the  sum¬ 
mer  months  and  .70  +.15  during  the  winter  months.  Although  the  secondary 
phase  correction  is  least  when  the  temperature  is  lowest  i winter;  anu 
corresponding  a  is  lowest,  the  variation  in  a  is  much  greater  in  the 
winter  than  in  the  summer.  For  this  reason,  the  temporal  variations  due 
to  changes  in  the  atmosphere  are  greatest  during  the  winter  months  in 
regions  which  experience  cold  temperatures. 

The  signal  tends  to  propagate  with  a  greater  average  velocity  in 
cold  weather.  Also,  during  cold  weather,  there  appears  to  be  a  correla¬ 
tion  between  loran  phase  and  temperature  that  does  not  exist  during 
periods  of  warmer  temperatures.  This  phenomenon  seems  to  be  related  to 
the  water  vapor  content  of  the  lower  atmosphere.  The  N-units  are  given 
by  equation  (4).  Although  the  first  term  (dry  term)  is  larger  than  the 
wet  term  (second  term),  the  change  of  N  due  to  the  water  vapor  alone, 
i.e.  3N/3e,  is  larger  than  3n/3t  and  3N/3P.  For  example,  for  a  temperature 
of  15°C,  a  pressure  of  1013  mb  near  ground  level,  and  a  relative  humidity  of 
60%  (i.e.,  e  ■  10  mb),  one  finds  =273  units  and  Nwefc  =45  units. 

However,  3N/ae  -  4.5/mb,  dN/ST-  -  1.26/  °K  and  3N/3P  *  0.27/mb.  Thus, 
as  the  temperature  decreases  far  below  the  freezing  temperature  of  water, 
the  partial  pressure  of  the  water  vapor  drops  significantly.  Thus,  it 
would  seem  that  as  long  as  the  temperature  is  above  the  freezing  point  of 
water,  the  wet  term  of  the  refractive  index  would  tend  to  keep  the  atmos¬ 
phere  homogeneous  over  the  lower  two  or  three  kilometers  of  altitude.  This 
region  is  considered  to  be  the  important  region  for  the  determination  of  the 
a  factor  which  in  turn  determines  the  precise  value  of  the  secondary  phase 
correction  (Reference  15).  During  the  periods  of  cold  temperatures,  there 
is  insufficient  humidity  to  keep  the  atmosphere  homogeneous  and  the  average 
value  of  a  assumes  a  value  of  a  J'.l.  Then,  during  extreme  cold,  the  a 
value  varies  from  a  s. 5  to  a  «ri.o. 
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4.2  Temporal  Variation  of  Propagation  Time  Due  to  Ground  Impedance 

Changes . 

Techniques  for  evaluating  the  ground  impedance  in  the  presence  of 
ground  horizons  at  various  depths  below  the  surface  have  been  given,  ref¬ 
erence  (6).  In  general,  the  impedance,  x,  is  a  complex  number, 

x  -  |x|  exp( j  Arg  x) , 

and  can  be  constructed  from  values  of  ground  conductivity,  o^,  c^,  03, 

...  ,  dielectric  constant,  c^f  Cj,  e^,  ...  ,  and  permeability,  ulf 
^2*  y3*  •••  »  *or  tf'e  various  ground  horizons,  provided  the 
electromagnetic  boundary  conditions  between  such  horizons  are  applied. 
Thus,  the  ground  over  which  the  electromagnetic  wave  propagates  is  in 
general  nonhomogeneous ,  not  only  in  the  horizontal  direction  but  also  in 
the  vertical  direction. 

The  electrical  ground  conductivity  or  resistivity  at  the  surface 
of  the  ground  (i.e.,  the  conductivity  of  the  surface  soil  horizon)  is  not 
always  the  only  important  parameter  in  the  determination  of  the  ground 
surface  impedance.  The  electromagnetic  wave  penetrates  the  ground  to  a 
considerable  depth  and  in  most  overland  situations  one  finds  that  the 
ground  wave  is  strongly  correlated  in  the  space  domain  with  subsurface 
horizons.  Thus,  the  geologic  age  and  characteristics  of  the  basement 
rock  influences  the  surface  impedance  and  hence  the  secondary  phase  cor¬ 
rection. 


4.2.1  Graphs  Depicting  Propagation  Time  Variation  Due  to  Ground 
Impedance  Changes. 


Figures  9,  10,  and  11  depart  the  effect  of  the  ground  surface 
impedance  upon  the  secondary  phase  correction.  The  values  of  effective 
ground  impedance  range  from 

x  »  .001055  exp (j. 7854),  (for  seawater) 
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SECONDARY  PHASE  CORRECTION,  MICROSECONDS 


4 


Figure  9  Secondary  Phase  Correction,  t  Microseconds  as  a  Function  of 
Distance,  d.  From  the  Transmitter  for  Various  Values  of  Ground 
Impedance  or  Conductivity,  a  *  0.65 
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Figure  10  Secondary  Phase  Correction,  t  Microseconds  as  a  Function  of 
Distance,  d,  From  the  Transmitter  for  Various  Values  of  Ground 
Impedance  or  Conductivity,  a  •  0.05 


SECONDARY  PHASE  CORRECTION,  MICROSECONDS 


to  an  extreme  value 


x  -  .08  exp ( j 1 .036) ,  (for  poor  earth). 

Vertical  lapse  factors  (a)  of  0.65,  0.85  and  1.20  are  used  respectively 
for  each  figure  (9,  10,  11).  These  figures  being  logarithmic  clearly 
show  the  effect  of  varying  impedance  at  short  distances  (<200  km). 

The  quantity,  |x|,  is  normalized  to  the  impedance  of  free  space 
(■  377  Ohms),  a  universal  constant  of  nature,  and  hence  is  actually 
dimensionless .  As  given  in  Figures  9,  10  and  11,  one  can  express  the 
effective  ground  impedance  in  terms  of  conductivity,  o,  mhos/m  and  dielectric 
constant  (relative  to  Cq  for  free  space,  e  *  The  perme¬ 

ability  vi,  is  of  little  significance  and  hence,  p  ■  po  »  4n  (10-7)  Henry/m. 
For  sea  water,  o  *  5  mhos/m.  In  the  extreme  case,  |x|  ■  .08,  and  the 
corresponding  conductivity  o  -  .00076.  Effective  surface  conductivities  of 
o  ■  .0001  have  not,  thus  far,  been  observed  in  the  natural  ground  of 
the  earth  at  frequencies  near  100  kHz.  Conceptually  it  is  difficult  to 
visualize  real  conditions  which  would  contribute  to  such  low  values  in 
effective  conductivity.  In  a  real  sense,  the  ground  surface  effective 
conductivity,  as  seen  by  the  electromagnetic  wave,  is  an  aggregate  of 
conductivities  contr ibuted  by  the  various  surface  and  subsurface  horizons 
extending  downward  to  an  order  of  the  skin  depth  for  the  wave.  The  con¬ 
tributing  biases  from  subsurface  rocks  and  soils  is  sufficient  so  as  to 
maintain  physical  boundary  conditions  at  the  higher  minimum  values. 
Physically,  it  is  highly  improbable  for  contributing  sources  to  permit 
such  low  values  of  effective  conductivity.  Accordingly,  values  of  0  - 
.0001  can  be  considered  nonphysical  in  so  far  as  the  propagation  of  radio 
waves  which  are  influenced  by  the  ground  is  concerned. 

The  impedance  concept  embraces  all  naturally  occurring  ground 
electrical  properties  which  have  been  observed  for  Loran-C  frequencies. 

Hence,  this  concept  provides  a  compact  generalization  of  physically 
realizable  electrical  properties.  Here  again  and  in  this  concept,  ef¬ 
fective  conductivities  like  a  -  .0001  are  nonphysical  and  have  no  mean¬ 
ing  for  Loran-C.  In  fact,  the  curves  shown  in  Figures  9,  10  and  11,  when 
plotted  parametric  in  0,  other  parameters  being  constant,  cross  over  as  the 
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the  value  o  -  .0001  is  approached.  Hence  the  secondary  phase  correc¬ 
tion  would  decrease  instead  of  increasing  as  a  function  of  the  corre¬ 
sponding  impedance.  To  avoid  this  nonphysical  domain  of  the  complex  im¬ 
pedance  plane  at  Loran-C  frequencies  it  is  necessary  to  move  the  phase 
angle  of  the  impedance,  Arg  x  -  it/4,  to  the  right  in  the  complex  im¬ 
pedance  plane.  This  results,  in  effect,  in  a  dielectric  constant  change 
as  shown  in  Figures  9,  10  and  11.  The  curves  on  these  figures  represent 
physically  realizable  values  of  ground  impedance  at  100  kHz  based  upon 
Loran-C  observations  in  the  continental  United  States,  Europe,  the  Medi- 
terrean  and  S.E.  Asia.  These  values  of  impedance  can  be  regarded  as  ef¬ 
fective  values  over  nonhomogeneous  and  irregular  ground.  Thus,  in  lieu 
of  the  use  of  the  propagation  simulation  for  nonhomogeneous  and  irregular 
ground.  Reference  (13) ,  practical  estimates  of  Loran-C  secondary  phase 
corrections  can  be  obtained  from  Figures  9,  10  and  11  and  experimental 
data  or  experience  with  the  Loran-C  system.  The  experimenter  and  naviga¬ 
tor  in  the  field  can  gain  through  experience  the  ability  to  estimate  the 
ground  impedance  for  any  particular  propagation  path.  Also,  the  physical 
realizability  of  a  Loran-C  observation  can  be  ascertained  immediately, 
especially  when  gross  errors  exist  in  the  data.  More  subtle  errors  can 
also  be  detected  with  experience.  Thus,  Figures  9,  10  and  11  present  a 
frame  work  of  theory  within  which  all  observations  are  subject  to  ex¬ 
planation.  However,  severely  anomalous  situations  will  require  the  in¬ 
tegral  equation/full  wave  for  homogeneous  and  irregular  ground.  Reference 
(13). 


4.2.2  Example  of  Temporal  Change  Due  to  Change  in  Ground  Impedance 

Consider  first  the  nature  of  the  ground  in  any  particular  loca¬ 
lity  as  it  changes  with  depth  below  the  surface  of  the  topmost  layer  of 
soil.  A  geoelectric  section  of  the  ground  between  the  surface  and  some 
depth  below  (on  the  order  of  skin  depth)  the  surface  of  the  topmost  soil 
horizon  should  be  considered  in  a  ground  surface  impedance  estimate. 

Thus,  the  waves  entering  the  ground  are  exponentially  attenuated  with 
depth,  or  distance  through  the  material  that  comprises  the  ground,  de¬ 
pending  upon  the  conductivity,  a,  or  the  resistivity,  p.  Although 
the  ground  is  layered  geologically,  such  layering  does  not  necessarily 
represent  an  electrical  boundary  like  much  clearer  air-ground  boundary  at 


the  surface  of  the  ground.  But  the  ground  does  vary  in  electrical  re 
sistivity,  pi#  dielectric  constant,  e^,  and  permeability,  i  -  1,  2,  3  . 
with  depth  below  the  surface.  Thus,  the  ground  is  anisotropic  in  the 
sense  that  the  average  resistivity,  is  different  in  a  direction 
parallel  to  the  surface  as  compared  to  the  vertical  direction  below  the 
surface.  The  basic  quantity  that  can  be  used  in  the  ground  wave  theory 
of  propagation  in  the  atmosphere  above  the  surface  of  the  ground  is  the 
ground  impedance,  x.  This  quantity  represents  a  value  of  impedance  re¬ 
flected  to  the  ground-air  surface  by  all  of  the  subsurface  material  of 
importance,  with  the  aid  of  this  impedance  an  electrical  boundary  condi¬ 
tion  for  the  ground  wave  can  be  established  at  the  physically  obvious 
air-ground  interface.  The  problem  is  how  to  model  the  region  below  this 
interface.  Experience  has  shown,  reference  (6)  that  a  three  layer  model 
is  usually  the  best  that  can  be  accomplished  in  most  land  areas  of  the 
world.  There  is  no  theoretical  limit  concerning  the  number  of  layers  or 
geoelectric  horizons  that  can  be  used.  The  limit  is  the  availability  of 
geophysical  data.  Thus,  in  the  three  layer  model  one  can  use  the  top 
soil,  the  subsoil  or  hardpan  (with  water  table)  and  the  geologic  struc¬ 
ture. 

This  simplified  model  can  be  used  to  illustrate  the  effect  of 
rainfall  upon  a  desert  propagation  path  between  Searchlight,  Nevada  and 
Ft.  Cronkhite,  California.  An  examination  of  the  geophysical  data,  re¬ 
ference  (6),  indicates  that  under  normal  dry  conditions  that  are  found  in 
the  desert,  the  topmost  soil  horizon  may  exhibit  great  resistivity  while 
the  second  horizon  is  saturated  by  the  water  table.  The  bedrock  is  often 
very  high  in  resistivi  y.  The  topmost  horizon  comprising  only  a  very  few 
meters  is  subject  to  sudden  changes  during  heavy  rainfall.  If  this  wea¬ 
ther  is  extensive  or  covers  a  large  portion  of  a  propagation  path,  chan¬ 
ges  due  to  ground  impedance  may  be  observable. 

Typical  conditions  along  the  propagation  path  over  the  desert  in 
terms  of  a  three  horizon  model  is  a  6  meter  top  horizon,  -  .0005,  below 
which  is  a  water  table  horizon,  02  .  .01,  with  a  thickness  of  10  meters. 

The  bedrock  can  be  assigned  a  conductivity,  -  .0005.  The  corresponding 
value  of  surface  impedance  for  desert  values  from  Table  4  iss 

x  ■  0.0325  exp^j.7646j  .  (case  A) 
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TABLE  4 

GROUND  IMPEDANCE  FOR  THREE  HORIZON  MODEL 
USING  TYPICAL  VALUES  FOR  DESERT  REGION 


i a 

n 

°2 

03 

W1 

W2 

|x| 

Arg  x 

■ 

.01 

.01 

6 

10 

.03320 

1.0230 

1 

.0005 

.01 

.0005 

6 

10 

.03254 

.7646 

.0005 

.01 

.0036 

6 

10 

.03312 

.9221 

.002 

.005 

.00054 

2 

10 

.0418 

.4939 

.01 

.01 

.04 

1 

2 

.0154 

.9596 

.067 

.067 

.067 

1.5 

10 

.00914 

.7847 

c 

.01 

.01 

.01 

1 

2 

.0236 

.7810 

0^  ■  Conductivity  of  top  most  horizon  of  soil 
a ^  •  Conductivity  of  second  soil  horizon 
Oj  ■  Conductivity  of  bedrock 

Wj  »  Thickness,  meters  of  top  most  horizon  of  soil 
W2  -  Thickness,  meters  of  second  soil  horizon 
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An  increase  in  the  conductivity  of  the  geologic  structure  to  a  value,  0 
■  .01  would  reflect  a  ground  surface  impedance  of  x  -  .0332  exp  j  1 . 023J 
(case  B) .  This  is  a  highly  inductive  ground  because  the  phase  angle  is 
considerably  greater  than  n/4  »  .7854.  Suppose  heavy  rainfall  causes 
the  topmost  horizon  to  reach  a  conductivity  o  ■  .01.  Then,  from  Table 
4  calculated  values  we  find: 

x  ■  . 023fi|exp  j.781oj  (case  C) 

Referencing  Figure  12,  using  Arg  x  -n/4  -  -.0044  radians  (case 
C)  *  tc  “  4.20  yis  (microseconds)  for  a  1000  km  propagation  path. 

Also,  x  ■  .033  using  Arg  x  -n/4  ■  .238  (case  B)  gives,  t  ■  5. 00  us 
or  a  change  of  about  0.80  us.  This  is  a  very  great  change  and  it  is 
clearly  physically  possible  but  perhaps  not  probable  under  normal  weather 
conditions. 

The  method  for  calculating  x  in  the  three  or  more  horizon  model 
is  given  in  reference  (6) .  It  is  interesting  to  note  that  a  path  in  the 
complex  impedance  plane  which  follows  |x|  ■  .001  and  Arg  x  -  n/4  - 
0  (or  slightly  negative)  up  to  a  value  of  |x|  »  .03  and  then  moves  to 
positive  values  Arg  x  -n/4  up  to  a  value  of  1.036  -  .7854  -  0.251 
radians  at  |x|  ■  .08  can  be  used  to  explain  most  Loran-C  phenomenon 
in  terms  of  secondary  phase  correction.  This,  of  course,  is  a  pragmatic 
simplification  and  leads  to  the  set  of  curves  given  in  Figures  9,  10  and 
11.  Thus,  Figures  9,  10  and  11  reduce  a  comparatively  complex  technology 
to  a  simplified  form  which  can  be  used  to  analyze  and  better  to  under¬ 
stand  Loran-C  coordinate  observations  or  Loran-C  ptopagation  time  obser¬ 
vations  (TOA's). 

Table  5  gives  values  of  t^,  where 

< 

»  3tc/3d,  ns/km 

parametric  in  the  impedance  magnitude,  |x|,  and  the  vertical  lapse 
factor  a.  The  corresponding  values  of  ground  effective  conductivity, 

0,  are  also  given.  These  values  are  derived  from  numerical  values  pre¬ 
sented  in  Figures  9,  10  and  11,  between  1000  and  1800  km. 
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TABLE  5 


RATE  OF  CHANGE  IN  SECONDARY  PHASE  CORRECTION,  t^.. 
NANOSECONDS  PER  KILOMETER,  FOR  VARIOUS  IMPEDANCE 
MAGNITUDES,  |x|  AND  VERTICAL  LAPSE  FACTORS,  a. 


1*1 

a 

1 

t^.ns/km 

a 

1*1  £ 

1 

t„, ns/km 

033 

.5 

3.497 

5 

.001055  .85 

2.233 

.55 

3.675 

.0557 

.01 

2.940 

.60 

3.844 

.0139 

.02 

3.701 

.65 

4.011 

.0051 

.033 

4.608 

.70 

4.168 

.0027 

.045 

5.420 

.75 

4.320 

.00076 

.08 

6.048 

.80 

4.466 

.85 

4.608 

.90 

4.746 

.95 

4.880 

1.00 

5.011 

It 

is  of  interest  to  note  that  a  ♦ 

18%  change  about  a 

-  .85, 

between  a 

■  .70  and  a 

-  1.00  produces  a  change  in  secondary  phase 

correction 

der ivative 

•  Atc  -  . 

843  ns/km. 

Also,  a  36%  change 

in 

the  impedance  |x|,  between  |x| 

*  .033  and 

.045  produces  a  change 

in  secondary  phase  correction 

derivative, 

Atc  -  .812  ns/km. 

The 

corresponding  change  in  the  ground  effective  conductivity  parameter,  0, 
is  47%.  it  thus  appears  that  the  impedance  magnitude  changes  and  the 
a-factor  changes  produce  the  same  order  of  magnitude  change  in  secon¬ 
dary  phase  correction  whilst  the  conductivity  must  change  about  11%  more 
to  produce  an  equivalent  change.  Thu 3,  a  is  less  sensitive  than 
| x |  and  a  in  producing  phase  variations.  In  this  we  have  assumed 
the  variations  are  about  average  values  a  -  .85  and/or  |x|  *  .033. 

Figures  9,  10  and  11  contain  information  to  that  available  on 
small  desk  top  calculator  computer  programs  using  polynomial 
representation.  These  computer  programs  can  be  used  to  derive  effective 
impedances  for  mixed  paths  using  the  Pressy,  Ashwell  Fowler  method, 
(Reference  24).  Table  6  gives  sample  calculations  for  three  segment 
propagation  paths  using  these  computer  techniques. 
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TABLE  6 

SAMPLE  CALCULATIONS  FOR  MIXED  PATHS 


Segment 

Segment 

Segment 

fcc 

No.  1 

No.  2 

No.  3 

Micro 

1 x | ,Ad(km) 

|x  ,Ad(km) 

|x  ,Ad(km) 

Seconds 

|x|  Effective 

.001,100 

.045,100 

.033,100 

1.525 

.0233 

.033,100 

.045,100 

.001,100 

1.525 

.0233 

.001,300 

.02,300 

.08,300 

5.126 

.0368 

.08,300 

.02,300 

.001,300 

5.126 

.0368 

Thus,  combining 

three  segments 

100km  in  length  each  yields  a 

secon- 

dary  phase  correction  of  1.525  microseconds  at  300 

km  or  an  effective 

impedance  of  .0233 

if  the  segment 

impedances  are  . 

001,  .045  and  . 

033, 

respectively.  Reversing  the  order,  of  coarse,  gives  the  same  result. 
Using  300  km  segments  and  values  of  each  segment  of  .001,  .02  and  .08 
respectively  gives  a  secondary  phase  correction  of  5.126  microseconds  and 
an  effective  impedance  value  of  .0368.  Thus,  all  of  the  mixed  path  wis¬ 
dom  including  the  net  path  impedance  values  are  contained  in  Figures  9, 

10  and  11  when  programmed  into  a  small  calculator. 
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5.  APPLICATIONS  OF  RESULTS 


5.1  Software  Program  for  User 

One  of  the  goals  of  thi3  development  is  a  portable  inexpensive 
machine  which  will  be  able  to  calculate  a  secondary  time  delay  as  fre¬ 
quently  as  the  situation  warrants  with  only  a  few  numbers  required  to  be 
supplied  by  the  user. 

The  general  calculation  for  the  secondary  phase  correction  of  the 
ground  wave  portion  of  a  IF  electromagnetic  signal  is  an  intricate  cal¬ 
culation  which  requires  machine  computation.  Many  experimentally  deter¬ 
mined  numerical  values  must  be  available  that  require  storage  allocation 
which  is  presently  only  available  on  large  machines.  By  using  appro¬ 
priate  average  values  for  some  of  the  required  data  the  calculation  can 
be  performed  on  a  smaller  machine  and  hopefully  result  in  correction  for 
the  secondary  time  delay  of  propagation  which  will  be  beneficial  as  a 
first  approximation  to  the  Loran-C  user.  A  FORTRAN- IV  program  has  been 
developed  which  performs  the  required  calculation  for  ranges  >  10Km 
with  a  50  nanosecond  accuracy  and  uses  only  a  moderate  amount  of  storage 
space.  A  flow  diagram  of  this  program  is  given  in  Figure  13.  The  Coast 
Guard  has  expressed  an  interest  in  having  such  a  program  for  use  on  the 
Hewlett-Packard  HP9825  desk-top  calculator.  A  listing  of  the  Software 
for  Functions  FAMP  and  TC  and  Subroutine  GROUNE  is  given  in  Appendix  A. 
This  listing  is  in  FORTRAN- IV  language  and  requires  translation  to 
Hewlett-Packard  language.  In  addition,  an  Executive  Main  Calling  program 
is  required. 

To  make  this  calculation  tractable,  parameters  in  the  general 
model,  which  introduces  effects  due  to  terrain  roughness,  soil  stratifi¬ 
cation,  soil  consistency,  geological  under layment  and  surface  elevation 
have  been  eliminated. 

The  parameters  of  importance  which  have  been  retained  are: 
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A)  Frequency  of  relevant  Loran-C  transmitter 

B)  Distance  of  separation  between  transmitter  and  user 

C)  User  elevation  above  surface. 

D)  Index  of  refraction  of  air  at  surface. 

E)  Impedance  of  surface. 

F)  Vertical  lapse  factor  for  index  of  refraction. 

The  first  three  parameters  are  scalar's.  Parameter  A  is  known 
exactly.  Parameters  B  and  C  can  usually  be  adaquately  determined  to  the 
first  order. 

The  remaining  three  parameters  are  of  a  more  complex  structure, 
being  non-uniform  spatially  and  temporally.  In  the  general  model  parame¬ 
ters  D  and  E,  which  have  a  two  dimensional  spatial  dependence,  are  de¬ 
fined  separately  at  each  node  of  a  specified  grid,  being  in  part  respon¬ 
sible  for  the  large  storage  requirements.  Parameter  F,  in  general,  has  a 
three  dimensional  spatial  variation  due  to  possible  inversion  layers  in 
the  atmosphere.  Numerical  data  describing  this  parameter,  since  it  is 
transient,  is  usually  not  available. 

Of  these  six  remaining  parameters  all  except  D  and  E  are  there¬ 
fore  the  same  as  in  the  general  model,  while  D  and  E  are  introduced  as 
single  scalar  quantities. 

The  index  of  refraction  of  the  air  at  the  surface  is  known  to  be 
a  function  of  temperature,  atmospheric  pressure  and  the  partial  pressure 
due  to  the  water  vapor  content  of  the  atmosphere.  The  empirically  ob¬ 
served  numerical  value  lies  between  1.0002  and  1.0004  and  in  the  abbre¬ 
viated  program,  the  parameter  has  been  built  in,  (therefore  need  not  be 
entered  by  the  user)  with  the  constant  value  of  1.000338.  This  is  the 
value  usually  used  for  the  average  propagation  time  calculation  of  the 
primary  wave. 

The  ground  impedance  values  corresponding  to  the  amplitude,  |x|, 
and  phase,  arg  x,  given  in  Table  1,  comprise  a  subtle  generalization  that 
permits  considerable  simplification  of  practical  Loran-C  predictions  and 
demonstrates  a  specialized  use  of  the  general  impedance  concept  to  Loran-C 
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that  has  been  found  to  be  useful  throughout  the  world.  The  application 
of  the  impedance  concept  is  based  on  observations  and  the  original 
observations  were  derived  in  References  17,  18,  19  and  20. This  concept 
has  been  successfully  used  in  a  variety  of  Loran-C  applications.  Since 
this  time,  further  use  of  the  impedance  technique  to  analyze  data  using 
small  programmable  computers  have  been  given  in  References  21,  22  and 
23.  The  method  is  readily  adapted  (as  demonstrated  in  the  above  refer¬ 
ences)  to  small,  portable  programmable  calculating  predictions. 

The  vertical  lapse  factor  is  known  to  vary  in  general  between 
0.60  and  1.20  due  to  weather  conditions.  For  the  abbreviated  model, 
values  between  0.70  and  1.0  are  acceptable  inputs. 

In  summary:  the  present  abbreviated  program  will  require  that 
the  user  enter  five  numbers  (A,  B,  D,  E,  F)  from  a  keyboard  entry  device, 
and  it  will  output  a  single  number,  the  secondary  phase  correction  in 
microseconds,  which  can  be  applied  as  a  first  approximation  to  correct 
for  position  location  on  the  usual  Loran-C  hyperbolic  grid  system. 


5.2  Assessment  of  Impact  on  Loran-C  Chain  Control 

'  Monitoring  and  fine-timing  adjustments  of  Loran-C  Signals  are 
performed  by  a  System  Area  Monitor  (SAM).  The  System  Area  Monitor  has 
receiving  equipment,  usually  located  remotely,  to  continuously  monitor 
the  relative  timing  between  master  and  secondary  station  transmissions. 
The  SAM's  normal  function  is  to  issue  timing  adjustments  to  the  secondary 
station  (Local  Phase  Adjustment,  or  LPA)  to  compensate  for  any  frequency 
offset  of  the  secondary  station  oscillator  with  respect  to  the  master 
station  oscillator,  and  includes  propagation  variations.  These  adjust¬ 
ments  maintain  the  phase  (cycle)  time-difference  at  the  assigned  value 
(Controlling  Standard  Time  Difference,  CSTD) . 

The  total  adjustment  thus  contains  both  the  oscillator  drift  and 
the  temporal  propagation  variations  along  the  path  between  the  secondary 
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stations  and  the  SAM.  This  approach  enhances  the  accuracy  of  a  user 
navigating  with  Loran-C  in  the  vicinity  of  SAM.  However,  because  the 
Loran-C  chains  and  SAM's  are  necessarily  land-based;  the  typical  user  is 
quite  far  away,  i.e.  on  the  ocean  or  in  harbors  or  estuaries.  Thus  the 
part  of  the  timing  adjustment  which  is  attributable  to  temporal  propaga¬ 
tion  variation  over  the  transmitter-SAM  path  is  not  applicable  over  the 
transmitter  to  user  psth  and  it3  inclusion  deteriorates  the  navigation 
accuracy  of  the  user. 

Therefore,  if  the  temporal  propagation  variations  at  the  SAM  can 
be  separated  from  the  total  timing  variation,  it  would  be  possible  to 
make  adjustments  which  would  compensate  for  oscillator  drift  alone. 
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6.  VALIDATION  OF  RESULTS 


In  order  to  validate  the  result*  presented  in  this  report,  it  in 
recommended  that  data  he  collected  at  a  site  (possibly  a  System  Area 
Monitor)  where  knowledge  of  the  electrical  character latics  of  the  prop¬ 
agation  path  Is  fairly  well  understood  and  in  a  geographic  location  where 
extreme  weather  condition?  are  experienced  (i.e..  Northeast  United 
States) . 


The  secondary  phase  correction  can  he  calculated  when  the  ground 
impedance,  x,  (or  conductivity  and  dielectric  constant)  of  the  propaga¬ 
tion  path  and  the  vertical  lapse  factor  of  the  refractive  index,  a,  are 
known.  If  the  electrical  character i st les  of  the  propagation  path  are 
known,  then  the  vertical  lapse  factor  becomes  the  most  important  param¬ 
eter  to  he  derived. 

TV  determine  a,  knowledge  of  the  temperature,  atmospheric  pres¬ 
sure  and  partial  pressure  of  water  vapor  are  needed  Iv^th  on  the  surface 
and  at  an  altitude  of  ~  2  km.  Optimally,  both  the  surface  and  altitude 
values  should  he  measured.  If  only  the  surface  values  are  available,  It 
will  be  necessary  to  develop  some  typical  vertical  profiles  for  various 
climactic  conditions.  Also,  the  immediate  history  of  the  surface  values 
of  T  and  e  may  be  used  to  determine  appropriate  lapse  rates. 

Using  these  measured  (or  dei ived)  values  of  a  along  with  con¬ 
tinuous  recordings  of  Loran  measurements  the  experimenter  can  relate  o 
with  measured  phase  changes  and  validate  the  results  presented  in  this 
report.  Ideally  all  measurements  should  I'e  made  throughout  the  J?4  hours 
of  the  day  and  throughout  the  four  seasons  of  the  year. 

Figure  14  illustrates  in  an  analytical  block  diagram  a  receiver 
processor  system  to  compensate  for  secondary  phase  error  due  to  meteoro¬ 
logical  effects  such  as  temperature  and  pressure. 

Although  the  block  diagram  shows  propagation  corrections  being 
sent  back  to  the  SAM  for  cesium  correct ions ,  this  system  could  very  well 
be  used  to  improve  the  fix  of  a  moving  platform. 


CESIUM 


Validation  Analytical  Block  Diagra* 


I'd  understand  how  this  system  operates,  consider  the  familiar 
expression  for  the  phase  of  a  low  frequency  radio  signal! 

♦  ■  7:  nad  ♦  ♦c  (Radians) 


This  expression  has  two  components,  #p  (primary  phase)  and  $c 
(secondary  phase).  4>p  la  calculated  from: 

♦p"  nj|  (Radians) 


and  4>c  ia  calculated  from: 


<k.a,  1/3c^3  t0  f 


(Radians) 


where  a  is  developed  from  a  meteorological  relationship  given  by 
Brenaner ,  H.,  1949  for  0°c  and  modified  by  Dr.  S.  N.  Samaddar  for 
general  uset 


a  *  1  ♦  a  x  10*<> 


where  e  •  partial  pressure  of  water  vapor  (millibars) 

T  •  Temperature  (°  Kelvin) 

P  ■  Total  pressure  in  millibars 
h  •  Height  in  meters 

The  primary  phase  (4>p)  can  be  calculated  in  a  straightforward 

way  at  a  fixed  SAM  site,  assuming  distance  (d)  is  known  accurately.  For 

a  moving  platform,  distance  becomes  a  dependent  variable  to  be 

calculated,  and  updated,  via  the  closed  loops  around  the  <J>_  and  $ 

P  c 

processes. 
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It  is  interesting  to  note  that  at  a  SAM,  e  and  T  may  not  be  known 
over  the  entire  path.  However,  dt/dh  and  de/dh  are  the  parameters  of 
importance  and  they  can  be  determined  at  the  SAM,  one  point  on  the  path. 
In  a  moving  vessel,  the  navigator  can  acquire  e  and  T  not  only  along  some 
segment  of  the  propagation  path,  but  also  as  a  function  of  time.  Imple¬ 
mentation  of  this  system  (Figure  14)  using  these  new  data  may  improve  his 
accuracy.  The  reader  is  directed  to  Section  4.1.2  for  a  more  detailed 
explanation  of  the  function  of  meteorological  effects  upon  Loran-C  accu¬ 
racy.  It  appears  that  predicted  accuracies  less  than  50  feet  may  be 
realisable  when  processors  sense  ambient  meteor logical  conditions  to 
fully  reap  the  benefits  of  Loran-C. 
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7.  CONCLUSIONS  AND  RECOMMENDATIONS 


7.1  Conclusions 


The  roost  important  parameter  of  concern  in  temporal  variations  is 
the  atmospheric  vertical  lapse  factor  ot.  This  factor  can  be  related  to 
the  surface  refractive  index  na,  which  directly  affects  the  primary 
wave  phase  correction.  Due  to  its  correlation  with  the  a  factor,  n 

3 

affects  the  secondary  phase  correction  only  indirectly.  Most  observable 
temporal  changes  due  to  the  atmosphere  can  be  taken  into  account  between 
values  of  a  -  0.60  to  1.20. 

The  changes  in  secondary  phase  correction  are  proportional  to 
A(cx^/3)  and  Ad  (increments  of  two  thirds  power  of  alpha  and  distance). 

The  effect  of  Aa2/3  increments  is  more  pronounced  at  greater  dis¬ 
tances  because  of  the  linear  increase  with  distance  as  a  multiplicative 
factor.  However,  the  variations  in  a  are  greater  in  the  winter  months 
than  in  the  summer.  Temporal  variations  in  observations  along  a  1000  km 
path  in  the  eastern  U.S.  can  be  explained  by  an  a  value  of  1.0  ±0.2 
during  the  summer  months  and  .70+. 15  during  the  winter  months. 

The  secondary  phase  correction  is  most  sensitive  to  the  elec¬ 
trical  constants  of  the  ground  impedance,  x.  Values  of  the  magnitude 
of  the  impedance  range  between  (x[  *  0.001055  for  sea  water  to  |x|  ■ 

.08  for  typical  areas  of  the  continental  United  States,  Europe  and  parts 
of  Asia.  The  latter  value  of  impedance  is  inductive  in  nature  and  cor¬ 
responds  to  a  conductivity,  o  *  .00076  mhos/m  with  an  accordingly  ad¬ 
justed  value  of  dielectric  constant.  Temporal  variations  can  also  occur 
due  to  a  change  in  the  impedance  of  the  ground  due  to  precipitation. 
Although,  as  described  in  the  example  described  in  Section  4.2.2,  a 
change  of  0.8  us  is  theoretically  physically  possible,  it  is  not  con¬ 
sidered  probable  under  normal  weather  conditions. 

A  flow  diagram  for  a  computer  program  to  calculate  the  secondary 
phase  correction  is  given  in  Section  4.3.  This  program  provides  an 
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accuracy  of  50  nanoaaconds  and  a  listing  is  available  in  FORTRAN  IV  lan¬ 
guage,  (Appendix  A) . 


7 . 2  Recommends  t i ons 

It  is  recommended  that  a  variations  be  established  by  using 
meteorological  data  for  different  altitudes  at  any  given  location. 

Loran-C  data  should  be  analyzed  using  this  value  of  a. 

It  is  recommended  that  the  aforementioned  program  be  translated 
to  Hewlett-Packard  language  for  use  on  the  HP9825  desk  calculators  nor¬ 
mally  used  by  the  Coast  Guard.  This  will  aid  both  in  evaluating  the  ef¬ 
fects  of  temporal  induced  errors  on  Loran-C  accuracy,  and  in  implementing 
chain  control  techniques. 

As  mentioned  in  the  objective,  area  monitors  are  used  to  control 
the  Loran-C  chains  but  the  propagation  temporal  changes  removed  at  the 
geographic  location  of  the  monitor  are  not  appropriate  for  some  other 
location.  The  use  of  the  above  mentioned  computer  program  might  be  use¬ 
ful  for  the  SAM  for  calculating  (out  of  the  timing  corrections)  the  sec¬ 
ondary  phase  correction  for  his  unique  location.  This  program  might  also 
be  useful  for  the  user  for  calculating  the  secondary  phase  correction 
unique  to  his  location. 

Lastly,  it  is  recommended  that  a  study  be  made  to  investigate  the 
possible  use  in  the  future  of  the  NAVSTAR  Global  Positioning  System  (GPS) 
for  Loran-C  chain  control  (in  lieu  of  applying  propagation  corrections  to 
the  present  control  techniques) .  GPS  provides  extremely  accurate  three- 
dimensional  position  data,  velocity  information  and  system  time  to  a 
user.  The  purpose  of  the  recommended  study  is  to  determine  whether  this 
type  of  information  will  improve  upon  the  current  Loran-C  chain  control 
provided  by  the  SAM.  The  present  system  corrects  for  errors  only  in  the 
geographical  area  of  the  System  Area  Monitor  while  all  other  areas  become 
degraded.  Using  GPS,  all  of  the  chain  area  would  be  treated  equally  and 
improvement  in  overall  accuracy  might  be  realizable. 


7-2 


REFERENCE 


1  "National  Plan  for  Navigation",  U.S.  Dept.  of  Transportation, 

November  1977 

2.  Eaton,  R.M.  Mortiner,  A.R.,  Gray,  D.H.  "Accurate  Chart  Latticing  for 
Loran-C",  Canadian  Hydrographic  Service 

3.  Illgen,  J.D.,  Nelson,  L.W.,  "Loran-C  Signal  Analysis:  Time  Varying 
Components  of  Loran-C  Equipment  and  Propagation  Induced  Errors", 
Private  Communication,  1978 

4.  Johler,  J.R.,  Kellar,  W.J.,  Walters,  L.C.,  "Phase  of  the  Low  Radio 
frequency  Ground  Wave",  NBS  Circular  573,  1956 

5.  Johler,  J.R.,  R.H.  Doherty,  "Temporal  Effects  of  Natural  Origin  in 
Loran-C  Ground  Wave  Propagation",  CRPL;  Report  79-1,  1  May  1979 

6.  Johler,  J.R.  "Geophysical  and  Geological  Data  Rase  Evaluation  for 
Loran-C  Ground  Wave  Propagation  Medium",  CRPL;  Report  78-9,  1978 

7.  Samadder ,  S.N.  "The  Theory  of  Loran-C  Wave  Propagation  —  a 
Review"  (U.S.  Coast  Guard,  G-DOE-4,  TP54,  Washington,  D.C.  20590; 
presented  at  1978  Helsinki  Conference  of  LTRSI;  to  he  published). 

8.  Doherty,  R.H.  and  J.R.  Johler  "Meteorological  Influences  on  Loran-C 
Ground  Wave  Propagation";  Jour,  of  Atm.  and  Terr.  Phys.  37,  117  to 
1124,  1979 

9.  Doherty,  R.H.  and  J.R.  Johler  "Analysis  of  Ground  Wave  Temporal 
Variations  in  the  Loran-C  Radio  Navigation  System,"  OT  Technical 
Memorandum  76-222,  (U.S.  Dept,  of  Commerce,  OT/ITS,  325  Broadway, 
Boulder,  Co.  80303),  1976 

10.  Dean,  W.N.  "Diurnal  Variations  in  Loran-C  Groundvave  Propagation* 
(Magnavox  Government  and  Industrial  Electronics  Company,  Fort  Wayne, 
Indiana),  1978 

11.  Bean,  B.R.  and  G.J.  Dutton,  Radio  Meteorology,  NBS  Monograph  92, 
(Suptd.  of  Doc.,  U.S.  Gov.  Print.  Off.,  Washington,  D.C.  20402), 

1966 

12.  Doherty,  R.H.  "Spatial  and  Temporal  Electrical  Propagation 
Measurement",  AGARD  Conference  Proceedings  No.  144,  1974. 

13.  Johler,  J.R.  "Loran-C  Ground  Wave  Secondary  Phase  Corrections 
Over  Nonhomogeneous  and  Irregular  Ground  Using  Transient  Signal 
Propagation  Techniques",  CRPLi  Report  78-12,  (U.S.  Coast  Guard 
Contract  DOT-CG-8429323A,  U.S.  Coast  Guard  Headquarters,  Washington, 
D.C.  20590;  Prepared  by:  Colorado  Research  and  Prediction  Laboratory 
Inc.)  1979 

14.  Nelson,  L.  and  E.  Feniger  (1978),  "Loran-C  Signal  Analysis:  Stability 
Experiment  Summary  Status  Report",  GE  78  TW-56  (General  Electrie- 
TEMPO,  Center  for  Advanced  Studies,  Santa  Barbara,  California). 

R-l 


15. 


Samaddar,  S.N.,  "Weather  Effect  on  Loran-C  Propagation,"  U.S.  Coast 
Guard  (to  be  published),  1979. 

16.  Craig,  R.A.,  I.  Katz,  R.B.  Montgomery  and  P.J.  Rubens te in, 
"Meteorology  of  the  Refract 've  Problem,"  Chap.  3  of  Radiation  Lab 
Series,  Vol.  13,  Propagation  of  Short  Radio  Waves,  McGraw  Hill, 
Edited  by  D.E.  Kerr,  1951. 

17.  Johler ,  J.R.,  "Loran  Radio  Navigation  Over  Irregular,  Inhomogeneous 
Ground  With  Effective  Ground  Impedance  Maps,"  OT/Trer  22  (  Suptd. 
of  Doc.,  U.S.  Gov.  Print.  Off.,  Washington,  D.C.  20402),  1971. 

18.  Fredericks,  R.J.  M.W.  Bird,  K. A.  Wheeler,  "Recent  Advances  in 
Accurate  Airborne  Computation  of  Secondary  Phase  for  Loran-C 
Radio-Navigation,  "  (Lear  Stegler,  Inc.,  Instrument  Division,  4141 
Eastern  Avenue,  S.E.,  Grand  Rapids,  Michigan  49508),  1975. 

19.  Johler,  J.R.,  D.C.  Hyovalti  and  W.B.  Jones,  "Impedance  Maps  to 
Predict  the  Effect  of  Irregular,  Inhomogeneous  Ground  of  Loran-C/D 
Radio  Navigation  Systems,"  OT  Report  73-6  (Suptd.  of  Doc.,  U.S. 

Gov.  Print.  Off.,  Washington,  D.C.  20402),  1973. 

20.  Fredericks,  R.J.,  "A  Simplified  Algorithm  for  Improving  the 
Absolute  Accuracy  of  Loran-C  Radio  Navigation,"  Journal  of  ION, 
Winter  1974-1975,  1974. 

21.  Doherty,  R.J.  and  J.R.  Johler,  "Interpretation  of  West  Coast 
Loran-C  Spatial  Errors  Using  Programmable  Calculator  Analysis 
Techniques,"  7th  Annual  Technical  Symposium  of  the  Wild  Goose 
Association,  4  Townsend  Road,  Acton,  Mass.  01720,  1978. 

22.  Johler,  J.R.  and  R.H.  Doherty,  "Significance  of  Groundwave 
Propagation  Anomalies  Observed  During  a  Loran-C  Chain  Validation," 
OT  Technical  Memorandum  76-22,  (U.S.  Dept,  of  Commerce,  OT/ITS,  325 
Broadway,  Boulder,  Colorado  80302),  1976. 

23.  Doherty,  R.H.  and  J.R.  Johler,  "Calibration  and  Evaluation  of  the 
Unattended  Ft.  Hood,  Texas  Loran  C/D  Chain,"  OT  Technical 
Memorandum  76-21R  (U.S.  Dept,  of  Commerce,  OT/ITS,  325  Broadway, 
Boulder,  Colorado  80302),  1976. 

24.  Pressy,  B.G.,  G.E.  Ashwell  and  C.S.  Fowler,  "The  Measurement  of  the 
Phase  Velocity  of  Ground-wave  Propagation  at  Low  Frequencies  Over  a 
Land  Path,"  Proc.  I EE  100,  pt.  Ill,  73-84,  1953. 


L 


AD-A076  211  ANALYTICAL  SYSTEMS  ENGINEERING  CORP  BURLINGTON  MA  F/6  20/14 

LORAN-C  SYSTEM  DYNAMIC  MODEL!  TEMPORAL  PROPAGATION  VARIATION  ST—ETC(U> 
JUL  79  L  M  CAMPBELL  *  R  H  DOHERTY  00T-CG-829358-A 

UNCLASSIFIED  ASECR-79-107  USCG-D57-79  NL 


END 

DATE 

FILMED 

7-80 

OTIC 


SUPPLEMENTAF 

\ 


INFORMATION 


f\  o7  & 


J 

Errata  Sheet 
for 

Loran-C  System  Dynamic  Model  -  Temporal  Propagation  Variation  Study 
Report  No.  DOT-CG-D57-79/NTIS  AD-A076214 

Rage  4-1:  Change  to  read  as  follows.  Vertical  indicators  in  margin 

denote  where  changes  have  been  made,  "and  the  phase  changes 
associated  with  the  surface  refractive  index  is  simply 

A  V  d, 
c 

which  quantity  is  in  microseconds  if  d  is  in  kilometers  and  c  in 
kilometers  per  microsecond.  Therefore,  a  change  from 

n  »  1.0002  to  1.0004, 

ol 

(the  maximum  na  variation)  over  a  propagation  path  with  a  length  of 
1000  km  would  yield  a  phase  change  of  0.67  microseconds.  Normally  the 
change  in  average  value  (n  -  1.000338)  over  a  day  or  during  a 
weather  change  would  be  only  40  N-units  (N-units  «*  (na  -  1)  106). 
Therefore,  the  maximum  diurnal  or  weather  change  effect  that  could  be 
attributed  to  changes  in  na  fQr  a  iooo  km  propagation  path  would  be 
0.133  microseconds,  or  133  nanoseconds.  This  amount  of  variation  may  be 
nearly  typical  of  the  variations  observed  over  a  1000  km  propagation  path 
during  warm  weather.  However,  in  cold  weather,  the  observed  variations 
are  considerably  greater  in  magnitude  than  these  values,  reference  (8). 

As  a  result  the  overall  winter  weather  variations  observed  over  an 
approximately  1000  km  propagation  path  have  been  attributed  to  changes  in 
the  vertical  lapse  factor.” 

Page  3-7:  After  Eq.  (5) ,  change  "h  is  in  meters"  to  read  "h  is  in  100 
meters". 

dN 

Page  6-3:  Change  Eq.  for  to  readf 

* 

dN  f776  x  1.268  77.6  (.  ,  9620e)  dT  77.6  /,  4810\  del 

dh  ,  -  ♦  ^2  Y  +  T  j  dh  +  T  l1  “  T  /  dh  J 

Page  6-3 i  After  Eq.  for  change  "h  ■  Height  in  meters"  to  read 
"h  ■  Height  in  100  meters." 
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